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We numerically examine the properties of a two-dimensional system of particles which have com- 
peting long range repulsive and short range attractive interactions as a function of density and 
temperature. For increasing density, there are well defined transitions between a low density clump 
phase, an intermediate stripe phase, an anticlump phase, and a high density uniform phase. To 
characterize the transitions between these phases we propose several measures which take into ac- 
count the different length scales in the system. For increasing temperature, we find an intermediate 
phase that is liquid-like on the short length scale of interparticle spacing but solid-like on the larger 
length scale of the clump, stripe, or anticlump pattern. This intermediate phase persists over the 
widest temperature range in the stripe phase when the local particle lattice within an individual 
stripe melts well below the temperature at which the entire stripe structure breaks down, and is 
characterized by intra-stripe diffusion of particles without inter-stripe diffusion. This is followed 
at higher temperatures by the onset of inter-stripe diffusion in an anisotropic diffusion phase, and 
then by breakup of the stripe structure. We identify the transitions between these regimes through 
diffusion, specific heat, and energy fluctuation measurements, and find that within the intra-stripe 
liquid regime, the excess entropy goes into disordering the particle arrangements within the stripe 
rather than affecting the stripe structure itself. The clump and anticlump phases also show multiple 
temperature-induced diffusive regimes which are not as pronounced as those of the stripe phase. 

PACS numbers: 61.20.Gy 
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There has been growing interest recently in under- 
standing how collections of particles with competing 
interactions order or disorder under different condi- 
tions. A general understanding of such behaviors would 
have a profound impact on our ability to control self- 
organization and to tailor the functionality in systems 
of this class [II-Qj such as collections of particles with 
completing interactions composed of a long-range repul- 
sion and a short range attraction [^-[ioj ■ Specific systems 
where this type of interaction may be realized include 
charge orderin g in cuprate superconductors [llMlSf , mag- 
netic systems [ij, [l^ , bubble, stripe and clump phases 
in two-dimensional electron systems [la . [l7|. a nd dense 
nuclear matter which forms pasta phases |18| . There 
are also a number of soft matter systems which can be 
effectively modeled as assemblies of particles with a fi- 
nite range repulsive interaction that competes with a 
shorter range attractive interaction |19l-l2ll|. Recent sim- 
ulations of colloidal systems with this type of interac- 
tion [20] generated many of the same types of phases 
observed in models with longer range repulsive interac- 
tions |9|] , indicating that many of the features found in 
systems with long range repulsion and short range at- 
traction will be relevant for other systems with compet- 
ing interactions as long the repulsive interaction is longer 
range than the attractive interaction. Clump, stripe, and 
bubble phases have also been shown to occur in systems 
where the particle-particle interactions are strictly repul- 
sive provided that there are two length scales present in 
the repulsive interaction 0, Q . 

For systems with long-range repulsive interactions and 
short range attractive interactions, the onset of the differ- 
ent clump, stripe, and bubble phases can be controlled by 



directly varying the relative strength of the repulsive and 
attractive interactions at fixed density, or by holding the 
interactions fixed and varying the density of the system. 
In the clump phase where groups of particles form indi- 
vidual islands, the clumps organize into a crystal pattern 
at larger length scales due to the long ran ge repu lsive in- 
teraction between adjacent clumps [a S BUS [201 • Other 
phases include stripe or labyrinth structures and anti- 
clump phases in which an ordered array of voids forms, 
along with a uniform phase that appears at very high 
densities. It is not known whether the transitions be- 
tween the different phases are well defined, whether there 
are additional subphases, or what is the nature of the 
disordering or melting transitions of the phases in the 
presence of thermal fluctuations. 

Our previous work on the long-range competing inter- 
action system employed fixed density samples in which 
the relative attractive strength of the potential was var- 
ied [5[ . We showed that for narrow regions of parameter 
space, the system forms a labyrinth or disordered stripe 
phase. Ordered stripes can be generated when the sys- 
tem is driven in the presence of quenched disorder; in 
this case, the stripes are aligned in the direction of the 
applied drive 5]. In the absence of quenched disorder, al- 
though the particles move under an applied uniform drive 
they do not experience the local shearing forces which are 
necessary to induce the plastic deformations required to 
align the stripe pattern. Simulations of a stripe-forming 
colloidal system with intermediate-range repulsive forces 
indicate that if a non-uniform drive is applied so as to 
create a shear in the system, then plastic deformation 
can arise and the stripes align in the direction of the 
shear [2^. For systems with long range repulsive in- 



teractions, it is generally difficult to reach a completely 
ordered ground state due to the inherent heterogeneity 
of the system. For example, if one patch of the system 
has a higher particle density than neighboring patches, 
it is difficult for an individual particle to escape from 
the overly dense patch due to the short range attractive 
forces, and it is also difficult for a particle to move into 
an underdense patch of the system due to the repulsive 
long range interactions [y, Q . Colloidal simulations with 
shorter range repulsive interactions also produced simi- 
lar results and showed the formation of partially ordered 
states ^. 

In this work we simulate the different types of states 
that arise as a function of particle density in a system 
with long range repulsive and short range attractive in- 
teractions. We introduce several new diagnostics to char- 
acterize these states and the transitions between them. 
We also examine how the states disorder as a function of 
temperature. The stripe states are particularly useful in 
understanding the thermal disordering process since dif- 
fusion along the stripes can be distinguished readily from 
diffusion between the stripes. We show that the stripe 
phase exhibits four temperature induced regimes, includ- 
ing a low temperature frozen phase, a phase where the 
particles diffuse along each stripe but where there is no 
particle transport between stripes, a regime where there 
is anisotropic diffusion but the stripe structure persists, 
and a high temperature phase where the stripe structure 
is destroyed. The onsets of these different phases produce 
signatures in the energy measures and energy fluctua- 
tions. We argue that since the stripe state can show an 
intra-stripe melting transition, the stripe structure may 
remain stable up to higher temperatures than the other 
possible structures, since the excess fluctuations in the 
stripe state can first melt the particles within each stripe 
rather than the stripe structure itself. We also show that 
similar melting behaviors occur in the clump and anti- 
clump states. 



following classical equation of motion: 



I. SIMULATION 



We consider a two-dimensional Langevin simulation of 
a system of size L^ x Ly with periodic boundary condi- 
tions in the x and y-directions. In order to potentially 
accommodate a triangular lattice, we take Ly = L and 
Lx = 1.097L. Within the system we place N interacting 
particles at a density p — N/{LxLy). In previous simula- 
tions we held the system size fixed and either varied the 
density by changing N or varied the relative strength of 
the attractive and repulsive particle interactions. Here 
we hold the particle-particle interaction strength fixed 
and fix A^ = 380 but vary L in order to alter the density. 
The initial particle positions are obtained using simulated 
annealing. The dynamics of a single particle i obeys the 
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The particle-particle interaction has the following form: 






-K.Ri. 



(2) 



where Rij — |Ri — Rj| is the distance between particles 
located at Rj and Rj. The first term is a Coulomb re- 
pulsion, meaning that for very short distances the inter- 
action is repulsive and the particles will not collapse to a 
point. The magnitude of the fluctuations experienced by 
the particles during the course of the simulation is such 
that Rij never drops below the minimum value Rij =0.1. 
Due to the long range nature of the Coulomb interaction, 
it is not possible to cut off the interaction range, so we 
use a Lekner real-space summation method to compute 
the long range term |22h24|. The second term is a short 
range attraction of Yukawa form with an inverse screen- 
ing length of K = 1.0 and with B = 2.0. The last term 
in Eq. 1 is the random fluctuation force from finite tem- 
perature which has the following temporal correlations: 
(F^it)} = and {FT'{t)Fj{t')) = 2rjkBTS,,Sit ~ t'). 



II. PHASES AS A FUNCTION OF DENSITY 
A. Clump Phases 

In Fig. l(a-j) we plot the particle positions after sim- 
ulated annealing in the low density regime p = 0.04 to 
p = 0.22. At the lowest density p = 0.04 in Fig. 1(a), a 
mixture of clumps appears where there are two to three 
particles in each clump along with a smaller number of 
individual particles. The number of particles per clump 
increases with increasing p so that at p = 0.10, for exam- 
ple, there are six particles per clump on average. The de- 
crease in the number of particles per clump with decreas- 
ing density that we find at the lowest densities suggests 
that for even lower density than we have considered, the 
limit of only one particle per clump would occur and a 
triangular lattice would form. The system failed to form 
an ordered ground state at the lowest densities we consid- 
ered even for a very slow simulated annealing procedure. 
This results from the very strong clumping tendency in- 
duced by the short range attraction combined with the 
relatively large inter-clump distances at these densities. 
Particles are strongly attracted into clumps, but the dif- 
ference in energy at larger scales between having a clump 
with n particles or n -I- 1 particles is too small to permit 
the particles that have already fallen into a clump to es- 
cape from the strong short-range attraction of the clump 
in order to form clumps of uniform size. The resulting 
polydispersity in clump size prevents the clumps from 
ordering into a larger-scale lattice. 
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FIG. 1: Particle locations in the entire sample after simulated annealing for different densities in the clump regime. The scale 
of the X and y axes changes in each panel, (a) p = 0.04, (b) p = 0.06, (c) p = 0.08, (d) p = 0.10, (e) p = 0.12, (f) p = 0.14, (g) 
p = 0.16, (h) p = 0.18, (i) p = 0.20, and (j) p = 0.22. The number of particles per clump increases as the density increases. At 
the lowest densities the overall clump structure is disordered; however, for p > 0.2 the clumps form a triangular lattice. 



Since the long-range Coulomb repulsion favors a tri- 
angular ordering of the clumps, the problem of trying to 
order a mixture of clumps of slightly different sizes has 
similarities to the filling of magnetic flux quanta at sub- 
matching fields on a triangular pinning array in a super- 
conductor [25.] . At matching fillings where the number 
of flux quanta is an integer multiple / of the number of 
holes, such as at / = 2, / = 3, or / = n, with n integer, 
a triangular flux lattice forms; however, at non-integer 
fillings, most of the configurations are highly degenerate 
and disordered states appear, such as at / = n -I- 1/2. 
At other fields such as / = n/3, where an ordered tiling 
of the triangular lattice is possible, the ground state has 
several degenerate orientations, so any ordered state that 
forms is very likely to contain grain boundaries. Analo- 
gously, in our system a possible partial ordering could be 
expected for a low density filling where 1/3 of the clumps 
capture n -f 1 particles and the remaining 2/3 of the 
clumps contain n particles. In superconducting network 
systems of the type considered in Ref. 2b] , the triangular 
ordering favored by the repulsive vortex- vortex interac- 
tions is strongly enhanced by the ordered substrate. In 
our system, the clumps move over a substrate-free con- 
tinuum, which adds even more degeneracy to the sys- 
tem and makes it very difficult to observe ordered states. 
Even though there is a difference of only one or two par- 
ticles in the number of particles in neighboring clumps, 
the polydispersity of clump sizes is largest at the low- 
est densities we consider, where the charge ratio between 
neighboring clumps can be as large as 1 : 2 or 1 : 3. 
It is known that large polydispersity in the particle in- 
teractions for two dimensional Sjrotems can lead to the 
formation of a disordered state [2a] . 

At higher densities where the clumps contain larger 
numbers of particles, the difference in the number of 



particles in neighboring clumps is still only one or two 
particles. As a result, the effective polydispersity drops 
significantly and the system crosses over to a limit where 
each clump has close to the same charge. This tendency 
for the clumps to have more triangular ordering as the 
clumps grow in size is clearly observable for p > 0.2 in 
Fig. l(i,j) and for 0.24 < p < 0.27 in Fig. 2(a-c), where 
the center of masses of the clumps form a triangular lat- 
tice. In Ref. [20|, for colloidal systems with shorter range 
interactions, a larger polydispersity in the clump sizes 
appeared; however, the clumps formed an ordered tri- 
angular state whenever all of the clumps contained six 
or more particles. In the short range, purely repulsively 
interacting soft-sphere system of Ref. [3|, triangular or- 
dering also appeared even for polydisperse clumps as long 
as the number of particles per clump was six or greater. 
Simulations of a mixed dimer and trimer clump state in 
a system with a two-step repulsive interaction produced 
disordered states very similar to the ones we observe [31. 
These results indicate that clump states with small num- 
bers of particles are generally disordered independent of 
the exact details of the particle-particle interaction po- 
tential, while clump states with larger numbers of parti- 
cles per clump tend to form triangularly ordered struc- 
tures. 

Figures 1 and 2 indicate that the clump states form 
for p < 0.28. As the clusters grow in size with increasing 
density, we find that there is a tendency for clusters with 
certain "magic" numbers of particles to appear more of- 
ten than clusters with non-magic particle numbers. In 
studies of isolated clusters of particles with the same in- 
teraction potential, the clusters with n — 14 particles 
were particularly stable [10|. A simple counting shows 
that clusters of size n = 14 appear the most often at 
the highest clump densities of 0.24 < p < 0.27. The 
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FIG. 2: Particle locations in the entire sample after simulated annealing for different densities in the clump and stripe regimes. 
The scale of the x and y axes is not fixed in all panels, (a) p = 0.24, (b) p = 0.26, (c) p = 0.27, (d) p = 0.28, (a) p = 0.30, (f) 
p = 0.32, (g) p = 0.34, (h) p = 0.36, (i) p = 0.38, and (j) p = 0.40. The clump phase appears for p < 0.28. For 0.20 < p< 0.28, 
the clumps form a triangular lattice. For 0.28 < p < 0.44 a stripe phase occurs. The stripe structures have some density 
modulations along their length for p = 0.28 in panel (d) and p = 0.30 in panel (f), but are more uniform for the higher 
densities. The stripe phases at 0.32 < p < 0.36 in panels (f), (g), and (h) contain three rows of particles per stripe. 



largest clusters of size n = 16 occur at p = 0.27, while at 
p = 0.28 we observe the onset of stripe phases. 

B. Stripe Phases 

In previous work on striped phases, the stripes 
formed disordered labyrinth patterns when a substrate 
of quenched disorder was present [5|. In the absence of 
quenched disorder, aligned stripes may form depending 
on the rate of thermal quenching. For rapid quench rates, 
different domains of stripe orientation can arise leading to 
a partially ordered stripe state; however, for sufficiently 
slow quench rates, the lower energy aligned stripe states 
can form. The orientation of the stripes is affected by 
the periodic boundary conditions and by the number of 
particles within each stripe. 

At the lower-density end of the stripe regime, we ob- 
serve stripes that are composed of three rows of parti- 
cles which themselves form a sub-triangular lattice struc- 
ture within the stripe, such as at p = 0.32 in Fig. 2(f). 
For densities p < 0.32, the particles within the stripes 
are partially disordered and form elongated quasi-clump 
states along the length of the stripe, creating pinch-like 
structures of the type seen in Fig. 2(d) and (e) at p = 0.28 
and p — 0.30. For the higher densities of 0.32 < p < 0.34 
in Figs. 2(f-g), each stripe is nearly consistently three 
particles wide. As the density increases, bulges in the 
stripe that are four particles wide begin to form in a 
staggered pattern, as shown in Fig. 2(h,i) for p ~ 0.38 
and p = 0.40. At p = 0.38 in Fig. 2(i), the stripes be- 
come partially disordered when a portion of the stripes 
buckle outward to form regions that are four particles 
wide. This process continues for p = 0.40 and p = 0.42 



in Fig. 2(j) and Fig. 3(a). For p = 0.44, Fig. 3(b) shows 
that a large portion of the stripes are now four parti- 
cles wide. Above p = 0.44, the system transitions into 
the anticlump phase of the type shown in Fig. 3(c). For 
different values of the screening length 1/k than we con- 
sider here, other types of stripe phases containing only 
two rows or as many as five rows of particles could be 
realized; however, the general evolution of phases from 
clump to stripe to anticlump as a function of increasing 
density should hold for any choice of parameters. 



C. Anticlump and Uniform Phases 

Figure 3(c) shows the anticlump phase, composed of a 
triangular lattice of voids, which first forms at p = 0.46. 
At the low-density onset of the anticlump phase, the an- 
ticlump state is anisotropic and two of the six neck re- 
gions that surround each void are thinner than the other 
four, as seen in Fig. 3(c). The anticlump phases for 
0.47 < p < 0.58, shown in Fig. 3(d-j), stih have the 
triangular void ordering; however, the neck regions sur- 
rounding each void are now all of equal width so that 
the anisotropy found at p = 0.46 is absent. The par- 
ticles outside of the voids attempt to form a triangular 
lattice segment that is roughly three particles wide. For 
p — 0.56 and 0.58, shown in Fig. 3(i) and Fig. 3(j), the 
width of this triangular lattice grows to be four particles 
wide before a transition to a slightly distorted uniform 
triangular lattice phase occurs at p = 0.6. This triangu- 
lar lattice uniform state persists for all p > 0.6 and the 
triangular ordering improves with increasing density. 
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FIG. 3: Particle positions in the entire sample after simulated annealing for different densities in the stripe and anticlump 
regimes. The scale of the x and y axes is nearly the same in all panels, (a) p = 0.42, (b) p = 0.44, (c) p = 0.46, (d) p = 0.47, 
(e) p = 0.48, (f) p = 0.50, (g) p = 0.52, (h) p = 0.54, (i) p = 0.56, and (j) p = 0.58. The stripe phase in panel (b) at p = 0.44 
contains stripes that are close to four particles wide. The onset of the anticlump phases appears in panel (c) at p = 0.46. At 
the low density side of the anticlump phase, the structure is anisotropic, as can be seen by the thinning of the necks between 
anticlumps along the x direction in panel (c). For densities p > 0.46, the anticlump phase is no longer anisotropic. 
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FIG. 5: (a) Average distance to closest neighbor (dmin) vs 
p. The jumps indicate transitions between the clump, stripe, 
anticlump, and uniform phases, (b) The variation a„ of the 
mean neighbor distance averaged over n neighbors for n = 3, 
4, 5, and 6 (from bottom to top). The dotted lines indicate 
the locations of the transitions seen in panel (a). Compared 
to (dmin), two additional features appear in (t„. There is a 
plateau in an for 0.2 < p < 0.28 corresponding to the regime 
where the clumps form a triangular lattice, while for p < 0.2 
the clumps are disordered. There is also a jump in a„ at 
p = 0.38 which corresponds to the density at which the stripe 
state changes from having stripes that are three particles wide 
to having stripes that are four particles wide. 



FIG. 4: Particle positions in the entire sample after simulated 
annealing for different densities in the uniform regime where 
the ordering is mostly triangular, (a) p — 0.60, (b) p — 0.64, 
(c) p = 0.68, and (d) p = 0.72. 



D. Phases and Neighbor Lengths 

To characterize the onset of the different phases, we 
compute the average closest neighbor distance (dmin) 
for the post-simulated annealing configuration at each 



density. To determine (dmin), we solve the all-nearest- 
neighbors problem with a simple algorithm. The dis- 
tance from each particle i to its closest neighbor at R^„ 

is dj,,i„ = |R, - TV^J; then, (dmin) = ^"^Ef f^min- In 
Fig. 5(a) wc plot (dmin) as a function of density p. For 
the lowest density clump phases, (dmin) is large due to 
the presence of many individual particles whose closest 
neighbor is located in a different clump. As p increases, 
(dmin) decreases as the particles within each clump are 
compressed closer together. The transition to the stripe 



phase at p = 0.28 is marked by a sudden jump in (dmin)- 
The stripes are able to pack the available space more 
efficiently than the clumps while maintaining a stripe- 
stripe spacing similar to the clump-clump spacing at 
the end of the clump phase. As a result, the particles 
within the stripes are able to decompress and shift to a 
larger average spacing as they form a partial triangular 
lattice within each stripe. The neighbor spacing (dmin) 
decreases with increasing p throughout the stripe phase 
as the particles within the stripe compress in order to 
maintain the stripe structure. The transition into the 
anticlump state ai p — 0.46 is again accompanied by a 
sudden increase of (dmin) as the particles expand in or- 
der to take advantage of the additional area available to 
them in the anticlump structure. Within the anticlump 
phase, (dmin) decreases with increasing p as the particles 
are compressed into the diminishing space between the 
voids. The transition into the uniform phase appears as 
a large jump in (dmin) ai p = 0.6. This jump is signif- 
icantly larger than the other jumps in (dmin) since the 
particles are now able to fill the entire system evenly. As 
p continues to increase within the uniform phase, (dmin) 
decreases as the lattice constant of the triangular lat- 
tice shrinks. The sharpness of the transitions in (dmin) 
suggest that these transitions may be first order in na- 
ture. In addition to computing the nearest neighbor dis- 
tance dmin for each particle, we can determine the dis- 
tance dj ^ , rij g {l...n\) between a particle and each of 
its n^ neighbors as identified by a Voronoi construction 
[231 ■ The distribution P(d"^) of d"' at a particular value 
of p is bidisperse for p < 0.6 due to the presence of two 
length scales: the short spacing between neighbors that 
are within the same clump, stripe, or on the same side 
of a void, and the longer spacing between neighbors that 
are in adjacent clumps or stripes or are on opposite sides 
of a void. 

We compute the variation in particle spacing through- 
out the system at different values of p using a measure 
(7„ constructed in the following way: For each particle, 
we determine the mean distance d^ from that particle to 
its n closest neighbors: d^ = n'^J^d^' for rij ranging 
over the n closest neighbors only. We then determine the 
mean value /i„ of d„ averaged over all particles, /i„ — 
N~^J2i '^n: and finally obtain the standard deviation 

an of this quantity, cr„ = y (A^ - 1)"^ I]f (^« ~ M«)^- 
In Fig. 5(b) we plot cr„ as a function of p for the num- 
ber of closest neighbors n = 3, 4, 5, and 6. At small 
values of p when the clumps contain only one or two par- 
ticles, the measure (t„ performs poorly for n > 2, but 
for larger clumps and in all the remaining phases, the 
measure is well defined for n < 7. The dashed lines in 
Fig. 5(b) indicate the locations of the transitions found 
in Fig. 5(a). As n increases, the magnitude of the vari- 
ance cr„ also increases for any given density p, which is 
expected since a larger number of neighbors are being 
averaged. Several of the features in (T„ correlate with the 
transitions found in (dmin) in Fig. 5(a); however, a„ cap- 
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FIG. 6: The cumulative particle displacements dx (dark lower 
line) and dy (light upper line) vs time at different tempera- 
tures for the stripe sample at p = 0.36 in which the stripes 
are aligned along the y direction, (a) At T = 0.031 both 
dx and dy are bounded, (b) At T = 0.78, the increasing dy 
indicates that diffusion is occurring in the y-direction along 
the length of the stripes, but the flat dx indicates that there 
is no diffusion transverse to the stripes, (c) At T = 2.0, the 
stripe structure is still present and the diffusion is anisotropic 
with dy > dx- There is now finite diffusion transverse to the 
stripes, (d) At T = 2.53, the stripe structure is melted and 
the diffusion is isotropic. 



tures at least two features which do not appear clearly 
in Fig. 5(a). The first is that for ti = 4, 5, and 6, an 
rapidly decreases with increasing p up to p = 0.2 and 
then reaches a constant plateau which persists until the 
transition into the stripe state at p = 0.28. The onset 
of the plateau at p = 0.2 is correlated with the den- 
sity at which the clumps start to form a triangular lat- 
tice, while for p < 0.2 the clump structure is much more 
disordered. This may indicate that the charge polydis- 
persity effect becomes small enough at p — 0.2 that a 
transition to a triangular clump lattice can occur. The 
second distinctive feature is the increase in an that oc- 
curs near p = 0.38 within the stripe regime. This feature 
is correlated with the density at which the stripes change 
from being an average of three particles wide to being an 
average of four particles wide. At the onset of the uni- 
form phase at p = 0.6, cr„ drops and collapses so that 
all (T„ take the same value even for different choices of n. 
The cr„ measure indicates that in addition to the primary 
clump-stripe-anticlump phases, there can be sub-phases 
associated with changes in the average stripe width or 
the onset of ordering of the clumps. 




FIG. 7: The particle positions (dots) and trajectories (lines) 
for a fixed period of time at different temperatures for the 
stripe phase at p = 0.36 highlighting the different diffusive 
regimes in Fig. 6. (a) At T = 0.031, the particles in the stripe 
are frozen, (b) At T — 0.78, the particles can diffuse along 
the stripes but there is no diffusion between stripes, (c) At 
T — 2.0, the stripe structure is still present and motion occurs 
both along the stripes and across the stripes. It remains easier 
for the particles to diffuse along the stripes than across the 
stripes, producing an anisotropic diffusion signature, (d) At 
T = 2.53, the stripe structure is destroyed and the particle 
motion is isotropic. 



III. TEMPERATURE INDUCED DISORDER 
A. Different Liquid Regimes in the Stripe Phase 

We first consider melting transitions in the stripe state 
at p = 0.36, where the stripes are aligned in the y- 
direction as shown in Fig. 2(h). To examine the effect 
of temperature, we calculate the particle displacements 
dx and dy as a function of time for this phase at different 
temperatures for both the x and y directions. 



N 



dx,y{t) = (1/A^) J2 \R^Ai) - R^Ato)\ 



(3) 



In an isotropic liquid, we expect d^ — dy. In Fig. 6 we 
plot dx and dy at temperatures T ~ 0.031, 0.78, 2.0, 
and 2.53. Figure 6(a) shows that at T = 0.031 the dis- 
placements are bounded at long times and that d^ and 
dy have approximately the same value. In Fig. 7(a), the 
particle trajectories during a fixed time at T = 0.031 in- 
dicate that the particles undergo small displacements due 
to the thermal fluctuations but that both the stripe or- 




FIG. 8: The value of the cumulative particle displacements dx 
(open squares) and dy (filled circles) after 2 x lO'' simulation 
time steps as a function of temperature T for the stripe system 
at p = 0.36. At low temperatures, there is a frozen regime 
with dx = dy\ at intermediate temperatures, an anisotropic 
liquid appears with dy > dx; and at high temperatures an 
isotropic liquid forms with dx = dy. 



dering and the crystalline ordering within the stripes are 
preserved, showing that the system is in a frozen state. 
We note that for much longer times than those plotted in 
Fig. 6(a), the values of dx,y do not increase but remain 
bounded at a nearly constant level. At T = 0.78, plotted 
in Fig. 6(b), dy grows without bound while dx saturates 
and is bounded. Fig. 7(b) shows that for T = 0.78 the 
particles readily diffuse along the stripes but there is no 
diffusion between the stripes. The fragmentary trian- 
gular particle ordering within the stripes is lost and an 
interesting state appears that acts like a liquid in one 
direction but like a solid in the other direction, which is 
characteristic of a smectic phase. For T = 2.0, Fig. 6(c) 
indicates that both dy and dx increase with time; how- 
ever, dy rises more rapidly indicating that the diffusion 
is anisotropic. In Fig. 7(c) the particle trajectories at 
T = 2.0 show that there is a large amount of motion 
along the stripes accompanied by a smaller number of 
places where the particles have hopped from one stripe to 
the next. The overall stripe structure remains intact. For 
T = 2.53, plotted in Fig. 6(d), dx and dy both increase at 
the same rate, indicating isotropic diffusion. Figure 7(d) 
shows that the stripe structure is destroyed at T = 2.53 
and the particle motion is uniform. Although the liq- 
uid state in Fig. 7(d) is isotropic, we observe short-lived 
clumps that do not occur in an isotropic liquid of purely 
repulsive particles |28j . 

To better characterize the onset of the different dif- 



fusive regimes we plot the value of dx and dy after 
2 X 10^ simulation time steps as a function of temper- 
ature T in Fig. 8. For T < 0.32, the particles within the 
stripe remain frozen and the diffusion is isotropic. For 
0.32 < T < 1.53, the particles can flow in a liquidlike 
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FIG. 9: (a) The standard deviation of the energy as vs T 
for the stripe phase at p = 0.36. The large peak at T = 2.4 
corresponds to the melting of the entire stripe structure. The 
smaller knee structure near T = 0.5 corresponds to the onset 
of melting of the particles within the stripes, (b) The specific 
heat C vs T for the same system shows a large peak at the 
stripe melting temperature and some features near the intra- 
stripe disordering transition. 



(a) 



(b) 



FIG. 10: The particle positions at one instant in time for 
the stripe phase at p = 0.36. The scale of the x and y axes 
changes in each panel, (a) T — 2.53 showing transient clump 
structures, (b) T = 6.13 showing that the transient clumping 
is greatly reduced. 



fashion along the length of the stripes but there is no 
diffusion transverse to the stripes. Within this phase, 
if we measure d^.y for increasingly long times, the value 
of dy continues to increase but dx remains at a plateau 
value corresponding to the width of the stripe, which is 
the maximum distance the particles can diffuse in the x 
direction in this regime. For 1.53 < T < 2.37, the stripe 
structure is still present although there is now some dif- 
fusion transverse to the stripes. The diffusion remains 
anisotropic with dy > dx- Finally for T > 2.37, the 
diffusion becomes isotropic when the stripe structure is 
destroyed. 

To further characterize the stripe melting, in Fig. 9(a) 
we plot ffE , the standard deviation of the energy fluctua- 
tions, as a function of temperature. We obtain this mea- 
surement by collecting a time series of the total particle- 
particle interaction energy E during 5 x 10^ simulation 
time steps while holding the system at a fixed temper- 
ature. We also compute the average energy E at each 
temperature and use this to construct the specific heat, 
C = dE/dT, plotted in Fig. 9(b). Both the fluctuations 



as and the specific heat C show a pronounced peak at 
T — 2.4 when the stripe structure melts into an isotropic 
liquid state. A knee structure appears in both measures 
near T = 0.3 at the onset of the intra-stripe liquid phase. 
The magnitude of the energy fluctuations as remains 
nearly constant for 0.3 < T < 1.4, inside the intra-stripe 
liquid. Both as and C gradually decrease with increasing 
temperature above the stripe melting temperature. For 
purely repulsive particle assemblies, the magnitude of the 
energy fluctuations falls off very rapidly above the melt- 
ing temperature. The slow decrease in as for T > 2.4 
in the liquid state of the stripe system reflects the pres- 
ence of short-lived clumps within the liquid which be- 
come smaller and even shorter-lived as T increases. In 
Fig. 10(a) we show a snapshot of the particle positions in 
the isotropic liquid phase at T = 2.53, showing that the 
particle positions remain highly inhomogeneous on short 
time scales and form clump-like structures that change 
rapidly with time. For contrast, in Fig. 10(b) we show 
a similar snapshot deeper into the isotropic liquid phase 
at T = 6.13, indicating that the transient clumping ef- 
fect is significantly reduced at higher temperatures. In 
Fig. 9(b), C has the same trends as as', however, the 
data is generally noisier. At the onset of the intra-stripe 
liquid near T — 0.3, a smaller peak-like structure appears 
in C, indicating an excess of entropy. 

An open question is whether the different regimes ob- 
served in the stripe state as a function of temperature 
are true phase transitions or only crossovers. The spe- 
cific heat data indicates that the melting transition into 
the fully isotropic liquid state is most likely a phase tran- 
sition associated with a change in symmetry. Simulations 
of similar stripe and pattern forming systems generally 
find that the transition from a stripe state to an isotropic 
state is first order [^, or in some cases is a dislocation 
unbinding transition [15J. Hysteresis measurements sug- 
gest that the transition in our system is first order; how- 
ever, our simulation size is not large enough to rule out a 
dislocation-mediated transition in the stripe phase. The 
much weaker peak in C at the onset to the intra-stripe 
liquid phase suggests that the onset of this regime may be 
a crossover. We note that in most of our realizations of 
stripe phases, there are usually some topological defects 
present in the triangular ordering of the particles within 
the stripe. These defects may permit low energy particle 
motions that could obscure a true peak in the specific 
heat at the onset of the intra-stripe liquid. In general, 
for a one-dimensional system there are no finite tempera- 
ture phase transitions; however, in the stripe system, the 
ordering within the stripe is two-dimensional in nature 
and the particles in neighboring stripes are coupled to 
one another. Thus, a truly frozen state could occur in our 
quasi-one-dimensional stripes at finite temperature if the 
particles within the stripes were completely ordered. The 
onset of the anisotropic diffusive regime is most likely a 
crossover since individual particles within a stripe can be 
thermally activated out of the stripe. Due to the struc- 
ture of the stripes, there should be a well-defined char- 




FIG. 11: The cumulative particle displacements D after 
2 X 10^ simulation time steps vs T for the clump phase at 
p = 0.26. The displacements increase rapidly once the clump 
structure is thermally destroyed. Unlike the stripe phase, 
the diffusion in the clump phase is isotropic. Also plotted is 
dD/dT vs T, showing a peak at the clump melting transition 
near T = 2.0 along with a small rounded peak near T = 0.25, 
which corresponds to the melting of the particles within each 
clump. 



acteristic energy required for these hops to occur. When 
the temperature corresponding to this energy is reached, 
the inter-stripe diffusion rapidly increases. 



B. Different Liquid Regimes Within the Clump 
and Anticlump Phases 

We next consider the same set of measurements for the 
clump phase aX p — 0.26. In Fig. 11 we plot D = d^ = dy, 
the isotropic particle displacements after 2 x 10^ simu- 
lation time steps, as a function of temperature in the 
clump phase. The diffusion is isotropic and Fig. 11 indi- 
cates that the displacements rapidly increase just below 
T — 2.0 with a cusp type feature appearing at T = 2.0. 
We also plot the derivative dD/dT, which passes through 
a peak just below T = 2.0 at the point when the clumps 
begin to break apart. There is also a small peak in 
dD/dT near T = 0.25 associated with the melting of 
the particles within each clump. If we perform the mea- 
surement over a longer period of time, the values of D 
for T < 1.125 remain unchanged while the values of D 
for T > 1.125 increase. 

In Fig. 12(a) we illustrate the particle positions and 
trajectories in the clump state a.t p — 0.26 for T — 0.031. 
Here the clump structure is ordered and the particles 
within each clump undergo little motion. At T — 0.78, 
shown in Fig. 12(b), the clump structure is still present 
but the particles within each clump are highly disordered 
and diffuse rapidly within the clumps. This intra-clump 
liquid state has many similarities to the intermediate 
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FIG. 12: The particle positions (dots) and trajectories (lines) 
for a fixed period of time at different temperatures for the 
clump phase at p = 0.26. (a) At T = 0.031 the clumps 
form a lattice and the particles within the clumps show little 
motion, (b) At T = 0.78, the clumps still form a lattice but 
the particles are diffusing within each clump, (c) At T = 1.28 
the clump structure persists but interclump diffusion is now 
occurring in addition to intraclump diffusion, (d) At T = 
3.125 the clump lattice structure is destroyed and diffusion 
occurs freely throughout the system. 



stripe liquid state. The overall structure of the system 
is ordered but subregions of the system are liquefied. In 
the stripe case, it was possible for the particles to diffuse 
along the full length of the stripes so that diffusion in this 
direction was unbounded. For the clump state, the dif- 
fusion in the intra-clump liquid always remains bounded 
by the finite clump size. In the stripe system the sharp- 
ness of the transition from the frozen stripe crystal to the 
intrastripe liquid is controlled by the amount of ordering 
within the stripe phase. In the clump system, a similar 
effect occurs. It was shown in previous work that there 
are certain magic sizes of individual clumps that produce 
highly ordered configurations of the particles within the 
clump [lO[. For example, a clump containing 14 particles 
can form a stable segment of a triangular lattice which 
is free of internal topological defects. Additionally, stud- 
ies of small numbers of particles in individual traps also 
indicate that the onset of rotational and translational 
disordering is affected by the number of particles in the 
trap and the degree of order in their arrangement |29l - [33 | . 
For our clump system, this implies that if all the clumps 
had the same number of particles and each clump was 
also well ordered internally, such as for n = 14 particles 
per clump, a well defined transition from a frozen clump 
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FIG. 13: (a) The standard deviation of the energy ge vs T 
for the clump phase at p = 0.26. There is a peak at T = 2.0 
at the clump melting transition, (b) The specific heat C vs 
T for the same system shows a similar peak at the clump 
melting temperature. In contrast to the stripe phase, there 
is no discernible feature in either quantity at the intra-clump 
melting transition. 



phase to an intra-clump liquid phase should occur. In our 
system, not all the clumps have the same number of par- 
ticles and some clumps contain particle numbers which 
produce less stable internal clump structures; therefore, 
our sample shows disorder at lower temperatures which 
smears out the transition to the intra-clump liquid. 

In Fig. 12(c) we plot the particle trajectories at T = 
1.28 in the clump state. The clump structures persist 
and the particles within each clump have liquefied; how- 
ever, there is now some diffusion of particles from clump 
to clump. This phase is similar to the stripe phase 
which showed anisotropic diffusion. In the clump sys- 
tem, the diffusion always remains isotropic. If we fol- 
low the trajectory of a single particle, we observe that 
it remains trapped within a single clump for a period of 
time, diffusing within that clump, before escaping from 
the clump and undergoing much more rapid diffusion un- 
til it is trapped by another clump. This type of mo- 
tion is very similar to the caging dynamics seen in cer- 
tain glassy systems where particles can locally diffuse in 
some region with an occasional long hop to another re- 
gion where the particle again undergoes localized diffu- 
sion. At very short times, the behavior appears purely 
diffusive when the particles have not had enough time 
to experience the confining boundaries of the clump. At 
intermediate times, the behavior is subdiffusive due to 
the clump confinement, while at very long times, the be- 
havior is diffusive since clump-clump diffusion can occur 
at this temperature. A similar behavior occurs for the 
stripe phases for diffusion in the direction perpendicular 
to the stripes. This behavior suggests that these regions 
may exhibit some glassy type behaviors, a topic which 
we will explore in a separate work. 

In Fig. 12(d) for T = 3.125, the clump lattice is de- 
stroyed and the system enters an isotropic liquid state. 
Just as in the stripe case, the liquid phase at tempera- 
tures not too far above the clump melting region shows 
short-lived density inhomogeneities which decrease in size 
and lifetime as the temperature increases. 

In Fig. 13(a) we plot aE versus T for the clump state 




FIG. 14: The cumulative particle displacements D after 2 x 
10^ simulation time steps vs T for the anticlump phase at 
p = 0.54. Here the onset of diffusion occurs at T = 0.28 
corresponding to the disordering of the particles within the 
non-void regions. This is followed by a second change in the 
behavior of D at T = 2.0, which corresponds to the melting 
or destruction of the void lattice. 




FIG. 15: The standard deviation of the energy ge vs T for 
the anticlump phase at p = 0.54. There is a peak at T = 2.0 
where the void lattice melts. There is also some indication of 
a change of slope near T = 0.28. (b) The specific heat C vs 
T for the same system shows a similar peak at T = 2.0. The 
specific heat falls off much more rapidly above the melting 
transition than in the clump or stripe phases due to the fact 
that the isotropic liquid for the anticlump phase does not 
possess the same type of density inhomogeneities found for 
the clump and stripe phases. 



at p = 0.26 and in Fig. 12(b) we plot the corresponding 
specific heat C vs T. Both quantities show a peak at T = 
2.0 where the clump structure breaks apart, followed by a 
slow decay at higher temperatures in the isotropic liquid 
regime. Neither quantity shows any noticeable features 
at the frozen clump solid to intra-clump liquid melting 
which is likely due to the fact that the clumps contain 
unequal numbers of particles. 

We next examine the effect of temperature on the an- 
ticlump phase at p = 0.54. In Fig. 14(a) we show the 
cumulative displacements D ~ dx = dy after 2 x 10^ 
time steps versus T. The displacements are isotropic and 
show a significant increase near T = 0.28 followed by a 
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FIG. 16: The particle positions (dots) and trajectories (lines) 
for a fixed period of time at different temperatures for the 
anticlump phase at p = 0.54. (a) At T = 0.125 the particles 
undergo little motion and the state is frozen, (b) At T = 0.72 
the void lattice persists but the particles have liquefied and 
diffuse throughout the non-void region, (c) At T = 1.71 the 
behavior is similar to that for T = 0.72 but with even more 
motion in the non-void regions, (d) At T = 2.42 the void lat- 
tice is completely destroyed and diffusion occurs throughout 
the sample. 



smaller change near T = 2.0. The triangular void crys- 
tal remains intact up to T = 2.0, but diffusion through 
the entire sample begins at significantly lower tempera- 
tures. In Fig. 15(a) we plot as for the anticlump phase 
at p = 0.54. There is a very sharp peak centered at 
T = 2.0 corresponding to the melting transition of the 
voids. There also is some indication of a change in slope 
near T = 0.28. In Fig. 15(b) we show the correspond- 
ing specific heat C versus T. There is a peak in C at 
T — 2.0 but there is no particular feature at T = 0.28. 
The specific heat falls off rapidly for T = 2.0 due to the 
fact that the liquid phase does not exhibit the same large 
short-lived density inhomogeneities found above melting 
in the stripe and clump phases. 

In Fig. 16 we illustrate the particle trajectories at dif- 
ferent temperatures in the anti-clump phase at p = 0.54. 
At T = 0.125, shown in Fig. 16(a), the particles undergo 
little motion indicating that the system is frozen. In 
Fig. 16(b) at T = 0.72, the triangular void structure re- 
mains present but there is now particle motion through 
the non-void regions, indicating that particles are able 
to diffuse through the entire system. In Fig. 16(c) at 
T ~ 1.71, the void structure persists but some particles 
can now move into the voids for short periods of time. 



(a) 



(b) 



FIG. 17: The particle positions at one instant in time at 
T = 2.42. (a) The anticlump phase at p = 0.54. (b) The 
clump phase at p = 0.26 



The behavior for 0.28 < T < 2.0 has the same features 
found in the intermediate liquid phase for the stripe state 
where the individual particles have liquefied but the over- 
all structure of the state remains solid. In Fig. 16(d) at 
T = 2.42, the void structure has melted and diffusion can 
occur anywhere in the system. 

In Fig. 17(a) we plot the instantaneous particle posi- 
tions in the anticlump phase for p = 0.54 at T = 2.42 
and in Fig. 17(b) we show the corresponding particle po- 
sitions at the same temperature in the clump phase at 
p = 0.26. Just as in the stripe phase, the clump phase 
shows strong density inhomogeneities in the liquid state 
that decrease in size and lifetime with increasing tem- 
perature. In the anticlump phase, the liquid state found 
above melting is much more homogeneous and as a re- 
sult the fiuctuations in the energy are much smaller for 
the anticlump system above melting than for the stripe 
or clump phases. We also find a uniform density liquid 
for the higher density uniform regime above its melting 
temperature (not shown). 

These results indicate that the clump, stripe, and an- 
ticlump phase all exhibit an intermediate fiuid phase, 
where the internal ordering of individual particles is lost 
and selective diffusion can occur while the overall struc- 
ture of the clump, stripe, or anticlump pattern remains 
intact. 



IV. DISCUSSION 

One of our main findings is that the competing long- 
range repulsive and short-range attractive interactions 
produce a system that can exhibit both solid and liq- 
uid behavior simultaneously. The liquid-like behavior 
occurs on the short length scale of the spacing between 
individual particles. Motion on this small length scale 
initiates at temperatures well below the temperature at 
which the large scale patterns formed by the particles 
break apart. Future studies of the impact of the co- 
existing liquid and solid behavior could include deter- 
mining how the existence of the two length scales affects 
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measurable time scales. For example, the short-scale liq- 
uid could be probed with short-time measurements while 
the larger-scale solid state could be accessed using longer 
time measurements. The liquid-like nature of the system 
on short length scales might also play a role in permit- 
ting the system to adapt itself easily to a substrate or 
to an applied external field. For example, the ordering 
of a large length scale might be disrupted by quenched 
disorder in a system that had no short length scale de- 
grees of freedom, whereas in our competing interaction 
system, the effect of the disorder could be compensated 
on the short length scale, permitting the large length 
scale ordering to be maintained intact. In essence, the 
existence of the short length scale degrees of freedom per- 
mits a screening of the disorder on the larger length scale. 
It would also be interesting to see which of the phases 
we observe, such as clump, stripe, or anticlump, is more 
adaptable to quenched disorder in this sense. Our sys- 
tem may have interesting connections to other types of 
disordered systems which exhibit liquid-like behaviors on 
some length and time scales but solid-like behaviors on 
other length and time scales. For example, glassy sys- 
tems have properties that are in some sense inverted from 
what we observe, since glasses have liquid-like behavior 
on long length and time scales but solid-like behavior on 
short time and length scales. The time and length scale 
separation produced by competing interactions has clear 
potential for controlling novel functionality in materials, 
particularly those materials which contain additional de- 
grees of freedom such as electronic, magnetic, supercon- 
ducting, optical, biological, and so forth. The competing 
interactions can achieve a low-dimensional confinement 
of active regions, such as our short length scale liquid, 
embedded in a self-consistent supporting matrix, such as 
our larger length scale solid. 



V. SUMMARY 

In summary, we have analyzed the different types of 
phases that can occur in a two-dimensional system with 
a competing long-range repulsion and a short range at- 
traction as a function of density and temperature. At 
zero temperature and low densities, the system forms a 
clump phase. The particles cluster together in clumps 
and, due to the long range repulsion, the clumps move 
as far apart from each other as possible and attempt to 
form a macroscopic triangular lattice. When there are 



few particles in each clump, the overall structure of the 
clump lattice is disordered due to the large charge poly- 
dispersity of the small clumps. As the number of parti- 
cles in each clump increases, this polydispersity effect is 
reduced and a triangular clump lattice appears. In ad- 
dition to the triangular ordering of the clumps, there is 
also a tendency for the particles within each clump to 
form ordered or partially ordered structures. For higher 
densities, there is a sharp transition into a stripe state. 
The particles within each stripe tend to have triangular 
or partial triangular ordering. At still higher densities, 
a transition into an anticlump phase containing a trian- 
gular lattice of voids occurs, and at even higher densities 
there is a transition into a uniform phase with partial tri- 
angular ordering. We introduced several new measures 
to characterize the onset of the different phases as a func- 
tion of density, including the calculation of the average 
distance to the closest neighbor. 

As a function of increasing temperature, we find that 
the stripe phase exhibits clear multiple-step melting. A 
frozen phase with little or no diffusion occurs at low 
temperature. With increasing temperature the system 
melts into an intra-stripe liquid, where particles diffuse 
along the length of the stripes but do not diffuse between 
stripes and the overall stripe ordering remains intact. At 
higher temperatures, there is an onset of stripe-to-stripe 
particle motion leading to anisotropic diffusion, and the 
stripe structure continues to persist. Finally, there is a 
transition to an isotropic liquid where the stripe structure 
is destroyed and the diffusion becomes isotropic. The on- 
set of these different regimes produces signatures in the 
specific heat and in the energy fluctuations. We argue 
that the existence of an intermediate intra-stripe liquid 
indicates that excess entropy can be put into the mo- 
tion of particles along the stripes rather than into de- 
struction of the overall stripe structure. The presence of 
excess entropy is manifested as features in the specific 
heat and energy fluctuations. We show that the clump 
and anticlump phases exhibit similar intermediate liq- 
uid regimes. In the clump phase, the particles remain 
confined within the clumps but are liquidlike and mobile 
inside the clumps, while the clump lattice remains solid. 
For the anticlump phase, the particles diffuse through- 
out the entire lattice while the void structure remains 
crystalline. 

This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. 
DE-AC52-06NA25396. 
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